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Abstract 

Quark mass effects in 0(a^) QCD corrections to the decay rates of inter- 
mediate Higgs bosons are studied. The total hadronic rate and the partial 
decay rate into bottom quarks are analyzed for the Standard (scalar) Higgs 
boson as well as for pseudoscalar Higgs bosons. The calculations of three 
different contributions are presented. First, the flavour singlet diagrams con- 
taining two closed fermion loops are calculated for a nonvanishing bottom 
mass in the heavy top limit. Their leading contribution, which is of the 
same order as the quasi-massless nonsinglet corrections, and the sublead- 
ing terms are found. Large logarithms arise due to the separation of the 
pure gluon final state from the bottom final states. Second, quadratic bot- 
tom mass corrections originating from nonsinglet diagrams are presented. 
Third, nonsinglet corrections induced by virtual heavy top quarks are cal- 
culated in leading and subleading orders. It is demonstrated that, in order 
the first contribution numerically dominates over the second and the 
third ones. 
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1 Introduction 



In the last years the Standard Model (SM) has faced a remarkable lot of experimental 
tests, which were performed at LEP and SLC with very high precision. Due to the 
overwhelming agreement between theory and experiments the description of the world 
of elementary particles through a SU(3)c x SU(2)l x U(1)y gauge theory has emerged 
as a profoundly tested and firmly established theoretical framework. Despite of these 
achievements the detailed nature of electroweak symmetry breaking is still waiting for 
experimental confirmation. The search for a physical Higgs boson and the study of its 
properties will be the main subject of future collider experiments. 

Standard Model properties of the Higgs boson have been discussed in many reviews 
(see for example jl], 0]). In the minimal SM one physical scalar Higgs boson is present 
as a remnant of the mechanism of mass generation. Particularly interesting for the 
observation of the Higgs boson with an intermediate mass Mjj < 2Myy is the dominant 
decay channel into a bottom pair H — > bb. The partial width T(H — >• bb) is significantly 
affected by QCD radiative corrections. First order 0(a s ) corrections including the full 
nib dependence were studied by several groups ||, |], ||, |{| [j], || . Second order corrections 
were calculated in the limit <C M%. Apart from the trivial overall factor due to 
the Yukawa coupling, corrections were obtained for otherwise massless quarks by ]|, 1C] 



and for a nonvanishing mass of the virtual top quark by [pi]] . Subleading quadratic mass 
corrections in the m^/M^- expansion were found in In this work we complete the 
discussion of quark mass effects in second order QCD corrections. Our calculations 
provide the so far missing contributions and in part also serve as a crosscheck for 
existing results. 

Theories beyond the SM are usually characterized by an enlarged Higgs sector and 
may allow for different quantum numbers of the Higgs bosons. For example, one of the 
most appealing extensions of the SM, the Minimal Supersymmetric Standard Model 
(MSSM), contains two complex isodoubletts with opposite hypercharge (see e.g. Q), re- 
sulting in five mass eigenstates of the Higgs fields: two scalar (CP-even) neutral H°, h°, 
one pseudoscalar (CP-odd) neutral A and two charged physical Higgs bosons. As 
a consequence QCD corrections to the fermionic decays of a pseudoscalar Higgs have 
been studied in the past in many works |14|, 0, [15|, . This was motivation enough 
to carry out our analysis of quark mass effects in decay rates of pseudoscalar Higgs 
bosons as well. Our formulae are tailored in such a way that they are immediately ap- 
plicable to the MSSM. However, the commitment towards this specific choice of model 
is limited by the fact that supersymmetric QCD with the exchange of gluinos is is not 
covered in this work. 

Let us start our analysis of QCD corrections to the decay rates of the neutral Higgs 
bosons into bottom quarks with the two point correlators 

U s / p (q) = i J dxe iqx (0\ T J s l p (x) J 5/Pt (0) |0) (1) 

of the scalar current J s = ^ j and the pseudoscalar current J p = fi'js'^ / for 
quarks with flavour / and mass mt, which are coupled to the scalar Higgs bosons (we 
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(b) 



Figure 1: Diagrams of the 0(a 2 ) singlet contribution. (Thick lines: top, thin lines: 
bottom, curly lines: gluon, dashed lines: Higgs) 

use the generic notation H for H° and h°) and the pseudoscalar one (A) respectively. 

The analysis of quark mass effects in second order QCD corrections is the main 
concern of this work. However, the calculation of 0(a 2 s ) corrections to the correlators 
keeping the exact quark mass dependence would be an enormous task. Fortunately only 
the first few terms of the expansion in the small parameter m 2 /s represent already a 
very good approximation for processes at high energies. The partial decay rates of the 
Higgs bosons into bottom quarks can therefore be written in the form 



(2) 



where all information is contained in the absorptive part of the corresponding current 
correlator: 
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The coefficients C 



S/P 

ff 



g S ^ P g S J P are the respective weights of the various cou- 



plings between Higgs bosons and fermions. For the MSSM they are given by g^p 



h° 



sin a I sin 0, gg° wn 

cot (5 and g^ own = tan/3 with a, (3 denoting the usual mixing angles (see e.g. BJ). All 



cos a I cos (3, g^p = cos a/ sin (3, g\ 



sin a I cos /3, g^. 
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(a) (b) 

Figure 2: Examples of 0(a>1) nonsinglet diagrams. 

formulae in this work are applicable for the Higgs decay in the mimimal SM with 
^ = 1,^ = 0. 

The outline of this work is as follows. The decay rate of the Higgs boson A involves 
traces over fermion loops which are separately coupled to a pseudoscalar current. For 
this case some technical details about the treatment of 75 in D = 4 — 2e dimensions are 
given in Section 2. Flavour singlet contributions to the partial decay rates are discussed 
in Section 3. The heavy triangle diagram of Figure la contains two closed fermion loops 
composed of a top quark and a bottom quark respectively. It gives rise to a contribution 
to Ar2, which is of the same order as the otherwise massless nonsinglet corrections. For 
the scalar Higgs boson it was presented in [17] some time ago. Subleading terms 



are also calculated in Section 3. The light triangle graph of Figure lb consisting only 
of bottom loops is suppressed by rnf/s. Its contribution to Afg , obtained by [12] for 



the scalar Higgs, is only part of the complete answer, since bottom mass corrections of 
the same order originate also from the heavy triangle diagram and have to be added. 
This is done for the scalar and the pseudoscalar case in Section 3. Again subleading 

(m) 

terms Al 2 are given. 

In Section 4 nonsinglet contributions (see Figure 2) to the partial Higgs decay rates 
are discussed. Quadratic mass corrections to R s l p in second order are presented in 
this section. Their contributions to Ar^ are obtained from a recent calculation of 
the corresponding current correlators [18| and agree with ||, H- Further more the 
calculation of the double bubble diagram of Figure 2a containing a virtual top quark 
loop is performed in the heavy top limit. For vanishing bottom mass the exact top 
dependence of this contribution was obtained in [|11[ . Our result for Ar2 agrees with 
the leading (power-suppressed) term in the Mjj/mf expansion of |Q. In addition we 
present the correction for non-vanishing bottom masses, thus completing the second 
order corrections Ar2 , AT 2 , , AT 2"^ • 
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In Section 5 the results obtained in the full six flavour theory are expressed in the 
framework of an effective theory with five active quark flavours. 

The numerical size of the corrections and conclusions are given in Section 6. 



2 The Treatment of 75 

The calculations of the two-point correlators are performed with dimensional regular- 
ization in the MS-scheme. For the pseudoscalar correlator we employ the definition of 
75 in D 7^ 4 dimensions as suggested by 't Hooft and Veltman |l9| . 

Following |^] we avoid an explicit separation of Lorentz indices into 4 and (D-4) 
dimensions and work, in D = 4 — 2e dimensions, with the generalized current 

where jl?^ = (7^7^ + j x Yl**Y - ^Yl^Y ~ 7 p 7 A 7 ! V)/4. Taking the limit 
D — » 4 at the end of the calculation, when the result is finite, the pseudoscalar current 
is recovered by 

J P 4^ W . (5) 
The corresponding generalized current correlator has the following form 



nXilte 2 ) = i / dxe iqx {<d\ T pW ( , )P [^VA 1(()) | Q) 

J (6) 

~ [p'^'X'] Ul{q ] 

where Jyj^rjL = (l/4!)det(g , Q , a /) with a = p,fi,v,X and a' = p',p,',u',\'. It is conve- 
nient to work with the contracted tensor nj^^J, which is related by 

U p = I (l + 8e + 52e 2 + — e 2 ) n[^ A ] °^ -H^] (7) 
to the pseudoscalar current correlator IT P . 



3 Singlet Contributions 

The flavour singlet diagrams are characterized by two closed quark loops which are 
connected by gluons only. The heavy triangle graph of Figure la containing a top 
as well as a bottom loop is calculated in the heavy top mass limit. In leading order 
it contributes to the quasi-massless corrections Ar2: The Yukawa couplings of the 
fermions yield a factor m^rrit- In addition each fermion trace produces a factor mj and 
mt respectively. Together with a power suppression of 1/m 2 on dimensional grounds 
all mass factors are combined to m 2 and are thus of the same order as the otherwise 
massless nonsinglet corrections. 
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This leading and the subleading power suppressed s/mf contributions are obtained 

, p4| , which effectively is an expansion 



by employing the hard mass procedure [21, 22 
in the inverse heavy mass l/m| and has sucessfully been used in previous works [25, 
|26|, for the calculation of singlet diagrams for the decay rate of the Z boson. All 



possible "hard" subgraphs containing the heavy particle are selected and expanded with 
respect to the external momenta and the small masses. This formal Taylor expansion 
is inserted as an effective vertex. 

Having integrated out the heavy top m 2 3> Mfj , m 2 the remaining diagram still 
contains two different scales Mfj 3> m?, where the Higgs mass comes into play through 
the external momentum of the propagator integral. In analogy to the hard mass pro- 
cedure the expansion in m 2 /q 2 is obtained through the hard momentum procedure, 
where all possible subgraphs, through which the heavy momentum q may be routed, 
are reduced to a dot and expanded with respect to m& and eventual small momenta as 
compared to q. In this way we have calculated the leading and next-to-leading order 
of the m\l ' s series. As a result of the expansions the computation is simplified due to 
a factorization of the integrals. The three loop diagram decomposes into the product 
of massive tadpole integrals and massless propagator integrals. The latter are com- 
puted with the help of the multiloop program MINCER 
manipulation language FORM p£ 



written in the symbolic 



For the light triangle diagram of Figure 2b containing two bottom loops we also 
perform an expansion in m|/s. As already stated in [12], the first nonvanishing contri- 



bution is suppressed by mf]/ s and reproduced by our calculation. It has to be combined 
with the corrections of the same order originating from the heavy triangle diagram. 
Additionally we find that the corresponding contribution from the pseudoscalar light 
triangle graph vanishes identically. 

The absorptive parts of the singlet diagrams (Figures la, lb) are given by 
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AT trian9le (A -» bb, gluons) 
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They are constituted by the sum of all possible cuts of the diagrams and represent the 
decays into two (bb,gg), three (bbg) and four particle (bbbb) final states. 

Note that the expressions (8) and (9) contain no explicit \i dependence. This is a 
direct consequence of the fact that the singlet diagrams contain no divergent subgraphs. 
This, in turn, means that their absorptive parts and, thus, (8), (9) have no dependence 
on the choice of the renormalization scheme. 

The contributions of additional "ultralight" singlet diagrams containing a bottom 
and another light (/ = u,d,s, c) quark loop deserve some extra discussion with respect 
to their assignment to partial decay rates into specific quark species. The absorptive 
parts of these diagrams comprise again two (bb, f f , gg), three (bbg,ffg) and four par- 
ticle (bbff) final states. The cuts with two fermion final states (bb(g), ff(g)) should be 
calculated separately and could uniquely be assigned to the partial rates T b i and 
respectively. The situation is different for the four fermion final state (bbff). An unam- 
biguous assignment to a specific partial rate is not possible in this case. The question 
to which partial rate this piece should most reasonably be counted must be decided 
according to the peculiarities of the experimental setup and identification methods. 
However, since the ultralight singlet contributions are proportional to m'jm^, they are 
much smaller than the already small vn\ contribution from the double bottom triangle 
graph. We therefore have adopted a pragmatic point of view and neglected the ultra- 
light terms in eq. (|8|). (For the decay of the pseudoscalar Higgs they vanish identically.) 
No complications of this kind arise for the total Higgs decay into hadrons. We agree to 
the result in (T^j for the light and ultralight singlet contributions to I^ad- However, in 
view of the above discussion we find that in this work the assignment of the triangle 
graphs to the partial rates T^j is not well motivated. The complete formula for IYad 
will be given in Section 5. 

In order to complete our calculation and to arrive at the singlet result for the decay 
rate T b i into bottom quarks one still needs to subtract from eqs. ® the pure gluon 
final state contribution. The decays H/A — * gg through a fermion loop can be found 



in the literature @ ||, @, ||, ||]. 
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Using the notation 
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The singlet contribution to the partial width of the Higgs bosons into bottom quarks 
is therefore given by 
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The distinction between bottom and pure gluon final states introduces large logarithms 
In 2 (ml/Mfj) and m\ In 4 (m 2 /M^) in the partial decay width r b g. Being due to the sub- 
tracted decay into two gluons, their physical origin might be traced back to kinematical 
configurations which correspond to collinear gluons with respect to the virtual bottom 
quark. This interesting feature deserves further analysis in the future. 



4 Nonsinglet Contributions 

Second order QCD corrections to the scalar and the pseudoscalar current correlators 
were calculated recently by one of the authors for the case when the external momentum 
is much larger than all relevant masses. The general case of nondiagonal currents with 
quarks of different masses was considered. A detailed description of this work will be 
published in a separate paper Jy| . The results can be applied to the Higgs decay rate 
in the special case of diagonal currents. 

Besides the three loop diagrams of Figure 2b, induced contribution from lower 
order nonsinglet graphs need to be taken into account, since they lead to second order 
corrections after renormalization of the coupling constant and the quark mass: 
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Z§ and Zfrjg are introduced 2Cj: 

8a 



7 p 



3 7T 



7T 



2 r 



1 



1 



18 27 



1 _ ^il 
7r e 



(16) 



+ 



7T 



1 /15 



1 \ 1/25 11 



12 1 ) e V 16 72 

Ordering the absorptive parts of the corresponding correlators in massless and 
sive contributions 
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where I = ln(M^ A //i 2 ). We find agreement with [12, 15| 



Note that in (|17|) and ( |18[ ) n'j stands for nf — 1. This is because only diagrams 
without virtual top lines are included there. It also implies that the renormalization of 
these contributions is to be done by means of (15-|l6|) with nf substituted by n'^. 

It should be stressed that we are still working in the full QCD including the top 
quark and nj = 6. The strong coupling constant and light quark masses are, thus, 
defined with respect to this theory. A transition to effective parameters relevant to the 
topless QCD with five active quark flavours will be performed in the next section. 

The last terms in eqs. (^9|),(|20|) arise from double bubble diagrams, where a light 
quark is running around a virtual fermion loop. Although for these diagrams a four 
fermion final state (bbff) is possible, the assignment to the partial decay rate into the 
primarily produced quark flavour T bb seems to be the natural choice. 

It remains to compute the contribution of the double bubble diagram Ar DB with 
a virtual top quark loop (see Figure 2a). Again the hard mass procedure can be used 
and leads to the result: 
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5 Transition to the Effective Parameters 



The results obtained above do not allow a naive decoupling of the top quark as they 
contain the log [i 2 /m 2 terms not suppressed by any inverse power of mj. From a practical 
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point of view such terms can hardly be considered as potentially large as fi is eventually 
set to be equal to the Higgs boson mass. Moreover, they may be summed up by using 
the methods of the effective field theory as it has been recently demonstrated on the 



example of top mass effects in the hadronic decay rate of the Z boson [34]. 

Still it seems to be reasonable to express all our results in terms of the effective 



and mf' appearing in the topless QCD with five quark flavours. The matching 
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The corresponding changes in our results can be summarized as follows: 



(24) 



As the singlet contributions start from order ft 2 , they do not change its functional 
form within our accuracy in a s . 

Nonsinglet contributions fli~9| ) and ( |20| ) are conveniently arranged not to change 
their functional form if new (ai 5 ^) 2 terms ( that are induced by terms of order 
(a^)° and after the replacements ( |l5| ) and ([U]) are done) are assigned to 
the diagrams with a virtual top loop. 



• The resulting expressions for the double bubble contributions read: 

3Gf 



AT DB (H -» bb) 



Ay/2- 

M 2 H 



mi 



+ 



TT 

"107 
675 

(5h2 
b ) 
2 

mi 



a.W 2 

TT 



(m, 



— In- 

45 

116 



M 



H 



mi 



+ 



8 , Ml 



In 



H 



(25) 



Ar DB (^4 -» bb) 
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The decoupling of the top quark holds true for leading as well as subleading 
bottom mass corrections. The leading (massless) order agrees with the expansion 
of 0. 
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6 Discussion 



In this section we combine the QCD corrections discussed above. The result is valid 
for intermediate Higgs masses (in particular Mjj < 2mj), since the singlet contribution 
was calculated in the heavy top limit. Everywhere below the notation a s = a^\fi), 
fhf = m^\jx) is understrood with / = u,d,s,c,b and fx = Mjj. Expressing the decay 
rate in the following form 
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the coefficients for the scalar Higgs partial width F(H —>■ bb) read 
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For the decay rate T(A — ► bb) of the pseudoscalar Higgs boson one has (/x 2 
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The formula for the decay rate of the standard scalar Higgs boson in the minimal SM 
may be obtained by setting C bb = C bb = 1. From eqs. (|28]) , (^9|) one derives 
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(33) 



For our numerical discussion we use as input parameters a top mass of mt = 176 
GeV and a bottom pole mass of m;, = 4.7 GeV. The latter translates into the running 
mass m b (Ml) = 2.84/2.75/2.69 GeV for Higgs masses of M H = 70/100/130 GeV. 
All other quarks are assumed to be massless. Based on Aqcd = 233 MeV one arrives 
at the following values for the strong coupling constant: a s (Ml) = 0.125/0.118/0.114 
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corresponding to the three differerent values of Mjj- The influence of the top quark 
induced contribution on the second order coefficient is significant. AI"2 is shifted to 
(36. 60/37. 39/38. \2)(a s /t:) 2 and is mainly due to the double triangle contribution with 
its mass logarithms. The quadratic bottom mass corrections prove to be rather small. 
They add a contribution of (—0.54/ — 0.32/ — 0.22) (a s / 'it) 2 to the second order result. 
Finally we reproduce the expression for the total hadronic width of the Higgs boson in 
the minimal SM: 
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(34) 

Contrary to the partial rate into bottom quarks the hadronic width is practically not 
affected by large logarithms ln(m 2 /Ml), since they vanish in the sum of bottom and 
gluon final states. Their absence is reflected in the numerical numbers for the second 
order coefficients. For the massless corrections they read (33. 55/33. 12/32. 82)/(a s /tt) 2 . 
The mass corrections are again small. They amount to (—0.19/— 0.085/— 0.048)(o^/7r) 2 . 
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